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Abstract 

In part 1 of this series compact integral representations had 
been obtained for the one-loop photon-graviton amplitude in- 
volving a charged spin or spin 1/2 particle in the loop and 
an arbitrary constant electromagnetic field. In this sequel, we 
study the structure and magnitude of the various polarization 
components of this amplitude on-shell. Explicit expressions are 
obtained for a number of limiting cases. 



1 Introduction 



As has been recognized many years ago [UEI the quantized Einstein- Maxwell 
theory predicts the process of photon-graviton conversion in an electromag- 
netic field. The tree level vertex for this amplitude is (see appendix A) 

^V(^/o + f" a F Ua ) ~ \^F^f af} . (1.1) 

Here hn V denotes the graviton, the photon, and the external field, k 
is the gravitational coupling constant. The corresponding photon-graviton 
vertex in momentum space reads — ^K,C^ u ' a , where 



(1.2) 



This interaction leads, assuming sufficient coherence of propagation, to photon- 
graviton oscillations which are analogous to the better-known neutrino flavour 
(3] and photon-axion oscillations jU EJ [6] . The true eigenstates of propaga- 
tion in a background field will in general be certain mixtures of photon and 
graviton states. Determining these eigenstates and their dispersion rela- 
tions requires, at tree level and in Fourier space, the diagonalization of the 



following matrix (see eqs. (A. 6), (A. 7), and (A. 8)): 



i_ K QWfi ^( r? M« r? "A +r?/i A 7? w_2 ?? ^ r? ^ +...) J \ h KX {k) ) 

Here ap{k) represents the photon and h K \{k) the graviton. For many cases 
of physical interest, this problem can be simplified assuming the field to be 
homogeneous or near-homogeneous, and the modified dispersion relations to 
be close to the vacuum ones. An efficient formalism for calculating the evo- 
lution of the photon-graviton or photon-axion system under these conditions 
was developed in [B] . 

Due to the smallness of the gravitational coupling k, the photon-graviton 
mixing case has received less attention than the photon-axion one. Never- 
theless, a number of authors have studied possible observable effects [2j El 
El |5J H2 [TOl [11]. The chances of observing this process in the laboratory 
appear very remote. As with other processes involving very small couplings, 



= 0. 
(1.3) 
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the natural setting is astrophysics where one can hope to at least partially 
compensate this smallness by large field strengths or exposure times. In 
principle, any process based on photon-axion conversion in a field (see, e.g., 
[T2"l [13] [151 US El EE] ) must have an analogue based on photon-graviton 
conversion. In [6] photon-graviton conversion was considered in a pulsar 
field, as well as in the galactic magnetic field, but the effect was found to 
be very small. Photon-graviton conversion in a primordial magnetic field 
has been proposed as a possible contribution to the cosmic microwave back- 
ground anisotropy [8j [9] . However, taking plasma effects into account ren- 
ders the effect negligible [10! . 

A natural enhancement of the photon-graviton oscillation occurs in the- 
ories with extra dimensions [19] due to the existence of an infinite tower 
of Kaluza-Klein gravitons. In [20] both the effect of the photon-graviton 
oscillation on the cosmic microwave background and the conversion in a 
pulsar background were reconsidered in this context, but the enhancement 
was found to be insufficient to lift these effects into the observable range. 

To our knowledge, the photon-graviton process has previously been stud- 
ied only at the tree level. In the first part of this series [21] (referred to as 
'part 1' in the following) we considered the one-loop corrections to this am- 
plitude due to massive charged spin (denoted n^") and spin 1/2 particles 
(denoted fl^'") in the loop (fig. 1). 




Figure 1: One-loop photon-graviton amplitude in a constant field. The 
double line represents the propagator of a charged scalar or spin \ particle 
in a constant field. 



Our motivation for considering this loop correction is twofold. First, 



the dispersion relation (1.3) will be modified in the following way by the 



one-loop contributions (see eqs. (A. 6), (A. 7), and (A. 8)): 



) 



= 0. (1.4) 

This equation involves the full one-loop photon-photon, photon-graviton 
and graviton-graviton amplitudes, computed in the constant external field, 
summed over all possible loop particles (including, e.g., photons and gluons 
for the graviton propagator). As is well-known, already in the pure QED 
case the one-loop corrected dispersion relation in an external field 



^k 2 -Wk v -W' v {k) a v {k) = (1.5) 



leads to highly nontrivial deviations from the vacuum case |22| 
EU [271 EH [291 [301 EH [321 E3]. We believe that it will be very instructive 



to generalize this study to the mixed photon-graviton system (1.4), partic- 
ularly considering the fact that in curved backgrounds superluminal phase 
velocities are known to occur [3U EH [36] . See also [37] for recent studies 
on the photon vacuum polarization in curved space at arbitrary frequencies 
and related investigations on a possible breakdown of microcausality. 



Second, while the tree level interaction term ( 1.1 ) depends linearly on the 
background field, the one-loop corrections depend nontrivially on the field 
strength as well as on the photon/graviton energy. Thus, although at linear 
order in F^iy the one-loop corrections are down by an explicit factor of o. 
compared to the tree level term, it is a priori conceivable that for sufficiently 
large fields and some range of photon energies the one-loop amplitudes would 
dominate over the tree level one. 

The plan of this paper is as follows. In section [2] we find choices of phys- 
ical polarizations well-adapted to the structure of the worldline parameter 
integrals obtained in part 1. In section [3] we specialize to the on-shell case, 
and introduce some convenient notation. In section[4]we consider the purely 
magnetic case, and explicitly evaluate the amplitude for three ranges of pa- 
rameters: in |4.1| we present a direct numerical evaluation for arbitrary field 



strength and photon/graviton energies below threshold; in 4.2 closed-form 



results are found for the zero energy limits; in 4.3 we consider the case of a 
weak field but arbitrary energies. 

It seems that photon-graviton conversion so far has been studied only 
for the magnetic field case. Although the physical relevance of this process 
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is even more hypothetical in the electric field case, in part 1 we kept the 
electric field component since, quite generally, in the worldline formalism 
calculations in a general electromagnetic field are not substantially more 
difficult than in a purely magnetic field [38, 39 . In this sequel, too, we 
shortly consider the electric field case in section [5] and present results for 
two cases where the electric result can easily be inferred from the magnetic 
one, namely the zero energy and weak field cases. We summarize our results 
in section |6j 



2 Polarization decomposition of the amplitude 

Let us now project the photon-graviton amplitude on physical polarizations. 
We need to choose two photon polarization vectors E\ i such that 



<^ = Vf = 0- (2.1) 

Similarly, we need two symmetric, traceless, and transverse graviton polar- 
ization tensors e^'- 



-i S — ZOij, £ i — £ i , £ ifl — U, K^ 1 £ i — U. 



(2.2) 



As we will see, substantial simplifications can be achieved by choosing polar- 
ization vectors which are adapted to the background field. In part 1 we had, 
following [38, 39 , written the worldline Green functions using the matrix 
basis Z±, 



^ v _ aF"" - bF^ > b F^ + aF^ 

+ a 2 + b 2 ' * a 2 + b 2 ' 

\ nv F^ X F X U - tfrf" f&V _ F^ X F{ + a 2 r]^ 



a 



2 



+ b 2 ' K -J a 2 + b 2 



(2.3) 

Here F^ v = ^s^apF^ is the dual field strength tensor 1 and a, b are related 
to the two standard Maxwell invariants by a 2 — b 2 = E> — E 2 , ab = E • B. 

1 We work in Minkowski space with rf v = diag( — h ++) and e 0123 = 1. 
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Note the following properties of this basis [381 EH] , 



-Z±, (2.4) 
Z + -Z_ = 0. (2.5) 

We now use this basis to define polarization vectors as follows |^] , 



el = 



(i± • kf 



Here the A±'s are normalization factors, 



A_ 



Z + ■ k) ■ (Z + ■ k) 



-k.{Z\)-k. 



-{Z- ■ k) ■ {Z- ■ k) = -iJk ■ (Z±) ■ k 



(2.6) 



(2.7) 



Together with the orthogonality relation (2.5 ) they ensure that £+,£- satisfy 



the conditions (2.1). Both A± can vanish; in this case the corresponding 



polarization vector is lightlike and cannot be normalized. We also note that 



Ai + A 2 



A: 



(2- 



(from (2.3)). The explicit form of e± in terms of E and B is rather compli- 
cated in the general case. However, as usual things simplify considerably if 
one specializes to a Lorentz system where E and B are both pointing along 
the positive z - axis, E = (0, 0, £/),B = (0,0,-B). This implies that a = B 
and b = E. Here the dependence of X±, e± on the field magnitudes drops 
out, leaving only a memory of the field direction: 



A+ = y(fci)2 + (A : 2)2^A ± , A_ = -i^(fe°) 2 - (k 3 ) 2 = A,, , 

(2.9) 

and 

2 For the photon polarizations this basis has been introduced in [32) with different 
conventions and notations (in particular, our b corresponds to — b there). 
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„ (O^-^O) = „ „ (fc 3 ,0,0,fc°) = „ 

(2.10) 

The subscripts _L, || refer to the field direction. Such a Lorentz system exists 
provided that E-B > 0. The case E-B < differs from this only by a parity 
transformation, but the case E • B = needs to be considered separately: 
This case can for E > B be transformed into the purely electric and for 
B > E into the purely magnetic field case. In both cases A± and e± are 



the same as in (2.9), (2.10) if the z-axis is chosen as the field direction. The 
remaining possibility is that E • B = and E = B, the "crossed field" case; 
here a = b = and the above basis cannot be used. 

To simplify further, without loss of generality we shall assume that the 
photon propagation is in the y — z plane (k 1 = 0,k 2 > 0). 

It will be useful to construct also the graviton polarizations using the 
same building blocks. We define 



) flV 



(2-11) 



Then the conditions (2.2) are fulfilled as a consequence of (2.1). 



This basis is extremely convenient since, when contracting lij^^p^ with 
these polarization vectors/tensors, many terms in the integral representa- 
tions obtained in part 1 (eq s. (3 .13) and (4.11) there) drop out on account of 
the orthogonality relation (2.5) and the antisymmetry of Z±. In particular, 
all terms involving a factor of Sb\2 • k or Sbi2 • k will vanish. 

We remark that a more standard, but less convenient, basis would be 
obtained by removing the time component of en by a longitudinal shift, 



?m = e v * * _ (o 12) 

II ii fcV(fc°) 2 - (£ 3 ) 2 ' 

The equivalence of these two choices may not seem obvious for the gravita- 
tional part, since neither C fJ,l, ' a nor Hs^'i^pin are transversal in the graviton 
indices. However, using the gravitational Ward identity ((A. 13) of part 1) 
it is easily shown that the shift makes no difference for the matrix elements. 
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In the following, we will denote 



C Aa = e^C» u ' a e a a (2.13) 

etc. where A = ©, <g> and a =_L, || . At the tree level, one obtains the following 
simple result: 



C® 1 - = -2BX ± , 

C® 11 = 2iEX ]] , 
C 0± = -2iEX h 

C®H = -2B\ ± . 

(2.14) 

The fact that for a purely magnetic (electric) field e® couples only to e 1 - 
(e") and e® only to e" (e -1 ) is a consequence of CP invariance (see [6]). It 
must therefore also hold for the loop corrections. 



For the one-loop correction, using the polarization choices (2.6), (2.11) 
in (3.13) resp. (4.11) of part 1 yields the following: 



frAa _ f^^-ism 2 ) z + z - 

ii scal — - . n / c 



n 



64-7T 2 Jo s } sinh(z + ) sinh(z_) 

x J*due- i °*J& 1 + ~ieC A * 

Aa _ eK 1^ ds ism 2 \ Z+Z-. 



spin 32tt 2 7 s [tanh(z+)tanh(z„) 

(2.15) 

where now z + = ieBs, z~ = —eEs and (repeated indices a,b, . . . are to be 
summed over ±) 



1 A a 

* = (2-16) 
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spin 



z +( j ^bii ^fh) z -(A B u A F11 ) (A^ l2 + A F22 ) X + 
+'iAg 12 \ + (S B ~ 12 S B1 2 — Sp 12 S F12 )\a 



^B12 



(A%12 + ^Fll) 2 (Ap 



F12) 



A 



+iA 



B12 



(^B12 + ^Fll)(^B12 + ^Fll) A Fl2 Ap 



F12 



A+A 2 , 



spin 



(S B12 S B12 Sp 12 S F12 )z- (A B12 + A F11 )z + (^A B12 + + z + A F12 A F12 



+iA B12 \-(S Fjl2 S B1 2 — S F12 S F12 )\a 



-%A 
-iA 



B12 



+ \2 



(A B12 + A F11 



AiA_ 



B12 



(A F12 

(^B12 + ^Fll)(^B12 + ^Fll) ~~ ^F12^F12 



AiA_ 



(2.18) 



The remaining components are obtained using the symmetry 



j: 



seal, spin 



j: 



seal, spin 



seal, spin V 
seal, spin V 1 



-), 



-)■ 



(2.19) 
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The integrands are written in terms of the standard worldline functions (see 
(3.21) and (4.9) of part 1) 



A ± 



A ± 



J B12 



A ± 

^B12 



A ± 

A B22 



A ± 

A B\2 



s 



± 

F12 



A ± 

n -F12 



A 



F22 



sinh(,z-i-(l — 2u)) 

sinh(2:-|-) 
cosh(z±(l -2u)) 1 

sinh(z±) z± 



coth(z-i-) 



1 



A ± - A ± 

A B\2 A B\\ 



cosh(z-|-(l — 2u)) — cosh(z-i-) 



sinh(z-i-) 



cosh(z-|-(l — 2u)) 

COsh(2:-|-) 

sinh(z-|-(l — 2u)) 

cosh(z-i-) 
tanh(z-i-) . 



(2.20) 



Note that, as usual in this formalism [38], 155] . the scalar loop integrands 



(2.17) are obtained from the spinor loop ones (2.18) simply by nullifying all 



fermionic worldline correlators S Fi j , A Fi - . 



3 On-shell amplitudes 

In vacuum at this point we would use the dispersion relation 



k 2 = 0. (3.1) 

The modifications of this relation due to gravitational corrections are not 
relevant for our present purposes, since they would produce terms of higher 
order in k. However, this is less clear for the field-induced corrections to 



the electromagnetic H fJ, ' u . The question of under which conditions (1.5) can 



still be well-approximated by (3.1) was, for the magnetic case, studied in 
|32j . There it was shown that this is the case at least for moderate fields 

2 

and frequencies, B < 0{B cr ) and u < 0{m), where B cr = — denotes the 
"critical" magnetic field strength (B cr = 4.4 x 10 9 T for electrons). Here 
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the restriction on the photon frequency is not very significant, since for fre- 
quencies beyond the pair creation threshold lo = 2m processes involving 
electron-positron pair creation become possible, and then are usually phys- 
ically more relevant than the dispersive processes which we are concerned 
with here. To the contrary, the bound on the field strength may pose a 
restriction for applications to magnetars which are believed to carry field 
strengths up to several orders of magnitude higher than B cr [40 . 

In the following, we will assume that the use of k 2 = is justi fied, and 
work out the consequences. Using k = lo = \k\ eqns. (2.9), (2. 10) simplify 
to (with k 1 = 0) 



e» = (0,1,0,0) 
1 



e 



ii 



sin ( 



cos 9, 0,0,1) 



(3.2) 



and 



Aj_ = lo sin 9 , 
Ay = — icusmO, 

(3.3) 

where 9 is the angle between the z-axis (the field direction) and the direction 
of the photon propagation. Since in this approximation the amplitude de- 
pends on lo and 9 only in the combination lo sin 9, there is no point in keeping 
the dependence on 9. We will therefore restrict ourselves in the following to 
the case sin 9 = 1 (propagation perpendicular to the field direction) . For a 
general field we are then left with the four parameters m, lo, B, and E. It 
will be convenient to work with the three dimensionless variables 



LO 

B 

E 



LO 

m ' 
eB 

m 2 
eE 
m 2 



B 

B cr 
E 
E 



(3.4) 
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Here E cr = — denotes the electric "critical" field strength (1.3 x W 18 V/m 
for electrons). Similarly, for the calculation of the integrals it will be useful 
to change to the dimensionless proper-time variable s = m 2 s. Moreover, as 
usual in this type of calculations we will change from u to v = 1 — 2u. 

Finally, since we wish to compare the one-loop and the tree level con- 
tributions, we normalize the former by the latter. Thus we will have to 
compute 

/ f\Aa \ 
jAa /- r> t?\ — T5„| scal.spin 



IL?« Wia (u,,B,E) EE ^[df^a) 

poo _^ pi 

= aRe — e~ lS / dv tt^ in (s , v , u; , B , E) 
Jo s Jo 

(3.5) 

with dimensionless integrands vr^aispin- 

4 The magnetic field case 

We specialize to the purely magnetic case, E = 0. Then z- = so that 



^B12 — u > 

^B12 = , 
^F12 = *-* ' 



(4.1) 



Using these identities in eqs.(2.17), (2.18), (2.19) one can immediately show 
that 



TTWII = tt»-L _ n 

seal, spin scal.spin ' 



(4.2) 



This is in accordance with the CP analysis mentioned above. 
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The integrands of the nonvanishing components are, using (3.2), (3.3), 



7T 



scal 



1 

8tt \ sinh(z^ 



exp 



A+ 1 

z+ 2 y 



2\\ 



2B 



7T. 



seal 



1 



8ir { sinh(^ 



exp 



A + 1 

z+ 2 V 



B12 



B12 



y 25 



+ A B12 v(S B12 - v t ^ 



(4.3) 



for the scalar case and 



7T, 



spin 



7T. 



spin 



1 / *+ 

4tt \ tanh(z H 



exp 



-4 + 



1 



B12 + ± (1 _ u 2 ) \^ 

^+ 2 V V 2j B 



(^B12) 2 (*5'f12) 2 + (^F12) 2 



(^B12 + ^Bll) (^B12 + ~ + ^JFll 



+^B12(('S'b12) 2 ( , S'b12) 2 (^-B12 + ^Bll) 2 + (^B12) 2 v ^BX2 + ^12) 



2 

25- 



' 3/ ' 

1 f ^4 



exp 



4 + 



B12 + ±( 1 _ t7 i J , 

2 V ^25 



47r [tanh(z_j 

J ^B12 — ^F12 ~~ ~(^B12 + ^Bll) + ^B12 ( V ^B12 ~~ ^F12 + ^ 



2 1 



25 



+ - 



(4.4) 



for the spinor case. Here z+ = iBs and S BV2 i ^b\2 e ^ c - are as m (2.20) with 
1 — 2u replaced by v. 



The resulting integrals (3.5) have a structure similar to the parameter 



integrals for the photon vacuum polarization in a constant magnetic field 
obtained by Tsai and Erber [29|. Evaluating this type of integral for the 
whole range of field strengths and photon energies is known to be difficult |29[ 
|31~1 W\\ W2\ 133] , and in fact appears to have never been done in full generality 
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for the photon-photon case. For energies below the pair creation energy 
a direct numerical calculation is unp r oblematic after rotating to Euclidean 



proper time, T = is. From eqs. (3. 5), (4. 3), (4. 4), (2. 20) it is evident that the 



v - integral is always finite, and the integrand of the resulting T - integral 
is exponentially falling for small Cj. The pair creation threshold u) cr can 
be recognized precisely through the fact that for Cj > u) cr the T - integral 
becomes divergent at large T. Thus, analyzing the asymptotic behaviour of 
the integrands for large T at fixed B one finds that, for the scalar case, the 
critical value Cj cr is the same for both polarization components, while for the 
spinor case it depends on the polarization: 



2^1 + B, 



Co m ■ = 2. 

cr,spm 



(4.5) 



This divergence signals the onset of real pair creation, although to study 
the pair creation process itself one would have to consider the full photon- 



graviton polarization matrix (see (1.4) above). All four threshold values (4.5) 



agree with what one finds for the corresponding photon-photon amplitudes 

[221 m\ M H2]. 

Beyond the critical energy Euclidean proper time cannot be used any 
more, while a numerical integration in the original proper - time variable s 
poses enormous difficulties, due to the combined effect of the oscillatory be- 
haviour of the universal exponential factor and the poles of the trigonometric 
pref actor functions. 

Thus we will restrict our investigation to the two regimes which have 
also been well-studied in the photon-photon case, namely (i) photon energies 
below the pair creation threshold with arbitrary B and (ii) arbitrary photon 
energies at low field strength. Special attention will be given to the limit of 
zero photon energy, since, as we will see, here one can obtain the amplitudes 
in closed form. 

4.1 Photon energies below threshold 

As m entio ned above, below the pair creation threshold the parameter inte- 



grals (3.5) are, after a Wick rotation T = is, suitable for a direct numerical 



evaluation at arbitrary magnetic field strength. Fig. [2] shows the results of 
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such an evaluation, using MATHEMATICA, for the spinor loop amplitude 
ratio f[® in (&, B). The amplitude ratio is shown for field strengths in the 
range B cr < B < 10 4 i? cr and photon/graviton energies Gj < 2. A global 
factor of has been omitted. 




Figure 2: Numerical plot of the amplitude n®"^((2>, B) for B between 1 and 
10 4 and uj < 2. A global factor of — r has been omitted. 

The corresponding plots for n^J n and n®^,n®[jj are very similiar for 
photon / graviton energies below the corresponding thresholds (the numerical 
integration becomes unstable for u too close to cj cr it). In particular, all four 
amplitude ratios display a logarithmic growth in B for large B 3> B cr . And 
indeed, for fixed uj < uj ct and B — » oo it is easy to show the following 
asymptotic behaviour of (the nonvanishing components of) in (d),5), 



■-spin 



a 

a 
3^ 



]n(B) 
ln(B) . 



(4.6) 



Thus this asymptotic behaviour is independent of the photon energy, as 
well as of the polarization choice. Moreover, it is easy to see that it relates 



directly to the renormalization terms in the integrands (4.3), (4. 4). Thus we 



recognize here a connection between the short-distance behaviour and the 
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strong-field limit of this amplitude which is familiar from the photon-photon 
well as other external field processes [43] |4*4"] . 



4.2 Zero photon energy at arbitrary magnetic field strength 

Closed- form results can be obtained in the zero energy limit. Setting Cj = 



in (4.3), (4. 4) one finds that the v integrals can be done analytically, yielding 



trigonometric functions of z + : 



dv n. 



dv 7T 



seal 
II 

seal 



57T 



3 o zlcoth z + — 62: + cothz + + z\ + 2 2 

— z + cosech z + 

2 22;_|_smhz + 3 

2 cosh z + 2 



dv 7T 



dv 7T, 



spin 

®IJ 
spin 



57T 

1 
47T 
1 

4ir 



sinh 2 z+ z + smh.z + 3. 

3 cosech 2 ,z + — coth z + ( 2z + cosech 2 z + -\ \ + 

V Z+ J 



2 cosech 



„ coth z+ 4 

2 + - 

z + 3 



(4.7) 



Thus the remaining proper-time integrals are of the same type as the stan- 
dard proper-time representations of the magnetic Euler-Heisenberg Lagrangians 
(see, e.g., [15]), 



^fcalC- 6 ) 



4 

m 


r°° ds 


16tt 2 


Jo s 3 


4 

m 


f°° ds 


8^ „ 





Bs (Bsf 
sin(Bs) 6 
Bs 



tan(Bs) 



(Scalar QED) 

(Spinor QED) 
(4.8) 



In fact, they can be expressed in terms of the derivatives of these Lagrangians 
as follows: 



nfeispinO^O^) 
ns eil,spin(^ = O,5) 



27ra/l d | d 2 

m 4 ^BdB dB 2 
Ana 1 d 

m 4 B dB 



C. 



EH 

seal, spin 



/"EH I r>\ 
^scaljSpinV / ■ 



(4.9) 
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The existence of this connection between the photon-graviton amplitudes 
and the Euler-Heisenberg Lagrangians is not an accident, 
been obtained before by Gies and Shaisultanov 



Eqs. (14.91) have 



using the linear coupling 
of the graviton to the Maxwell stress tensor (see eq. ( A.2| )), and exploiting 
the relation between the vacuum expectation value of the Maxwell tensor 
and the effective action, 



2 5T 



(4.10) 



Applying (4.9) to the explicit representations of the Euler-Heisenberg La- 



grangians in terms of the Hurwitz £ - function \A7\ |4"8| |4"9| |4"5] 



4vr 2 



1 m / 



2 2eB 
3 / m 2 \21 



1 m 



+ Ch -1,- + 



2 2e£ 



1 +ln( 

\2eBJ 



■'spin 



2vr 2 



A\2eB 

,2 



f , / m 2 \ / m z \ / m z \ 1 1/m \2l 

ln (2^)-l2 + 4(2^) } 



(4.11) 



(the prime on C,h refers to a derivative in the first variable) one obtains 



n s e p f n (o,5) 



nfi(o,B) 



8^ 



a 

87 



1 lul^.li 



-1,- + — - + 

2 2B' ^2 



6 , „/l 1 
5 v 2 25 



-iln(2^) - 2 ln(2m - i- 
B 3 B 2 
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l,'+i)-ilnr(UiWiln(2. 



2 25' 5 ^2 2B J B 



-]n(2B)- + - 
3 v ; 5 2 3 



47T 



iinr 

5 v 25 



£ 2 v 25' 



1 / 4vr 3 x 4, , A , 2 1 

ta ( ) + iM2B) + 



47T 



+^(2«)-i 



ilnrfi) - ihirB/Tr) 
5 ^2B J B 

2' 
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Here ifi(x) = T'(x)/T(x) is the digamma function. 



(4.12) 



4.3 Weak magnetic field and arbitrary photon energy 



We proceed to the case of arbitrary photon energy but weak magnetic field, 
i.e., B <C B CT . Our treatment of this case parallels the one introduced by 
Tsai and Erber for the photon- photon case [29] (see also [32 J. 

Assuming B <C 1, we can expand the trigonometric functions appearing 



in the integrands (4.3), (4.4), in the common exponential factor (which is the 

]) as well as in the prefactor functions, 



same as in the photon-photon case 
keeping only the terms of lowest order in z+. However, since we do not wish 

* 2 

to make any assumption on the ratio this truncation has to be done 

* 2 

separately for the terms with and without a factor of For example, the 



prefactor of 7r® al involves (see eq. (4.3)) 



(S 



B12> 



A 



B12 

2 



2 



1 . 



+ — ) + A+ 



B12 



+ \ ^ 2 



(^512) (^B12) V ^B12J 



+ 



tjj 



B 



48 



(v 2 -i) 2 (3- w 2 )4 + o(4) 



2B 
(4.13) 



Although the first term in the square bracket is O(zl), at this stage it cannot 
be neglected over the leading \ + \ term. Similarly, the expansion of the 
universal exponent in eqs.(4.3), (4.4) yields 



* + (^ + i(l-» 2 ))^ = i(l-»V|4+0(4). (4.14) 

Here the leading order term, although of order O(z^), cannot be neglected 
with respect to the exponent of the global factor e~ ls . 

After performing these truncations, the amplitude ratios turn out to 
depend on B and u only in the combination BCj. This motivates the intro- 
duction of a new parameter A, 



-BCj. 



(4.15) 



Moreover, for the following it will be useful to interchange the orders of 
integrations in (3.5 ), and perform a v - dependent change of variables of the 
global proper-time variable from s to x, 
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where 



s = px 



(4.16) 



( 6 y 



(4.17) 



After this truncation and change of variables, the amplitude ratios take the 
following form: 



* (u 2 + l) (e" iH - e"^) - 

flfUX) = -^Ref'dvf 00 - 
spm 4vr 7o 7o x 

^ (t> 2 - 3) (e- iS - e~ ipx ) + |(3 + t; 2 )x 3 e- iS 



e -^_ e -i^ _ ,3 _ v 2^ x 3 e -i= 



i-v 2 x 3 e lS 



a 



nr P in(A) = —Re / dv 



1 roc 



4tt J Jo x 



dx 



X - (v 2 - 3) (e~ iE - e~ ipx ) +4(3- v 2 )x 3 e~ iE 



(4.18) 



The common exponential factor e ia is now of the Airy form, 



px+ -x 



(4.19) 



Note that we have also rewritten the counterterms in a way which will be 
convenient in the following. 
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Now, the two different x - integrals appearing here can both be expressed 
in terms of the modified Airy function Gi(p), defined by 



If 00 / 1 

Gi(p) = — dxs'mipx + -; 

vr Jq V 3 



(4.20) 



Namely, one has 



flip) = 



r°°dx[ / i 3 \ , N - 

/ — cos [px H — x ) — cos(px) 
Jo x L V 3 J 

dx \ ( 1 3\ f i i ■ ( i 1 ' 

— cos -x — x / a/3 sin o x H — x' 
x L v 3 / Jo v 3 

vr / P dp' Gi(p') + ln / 9 + -7 + -ln 3, 
JO 3 3 



cos(px) 



(4.21) 



and 



/ 2 (p) 



dx x sin 



(px + * x 3 ) = -vr G*"(p) = 1 - vrp Gi(p) . 



(4.22) 



This brings us to our final result, 



n s e ii(A) 

n^in(A) 
nfp'L(A) 



5vr Jo 



5VT Jo 



a 
4vr 



a 
4vr 



dv 



1 -^ + l)f 1 (p)- 1 -(3-v 2 )f 2 (p) 

l -(y' + l)h{p)- 2 -V 2 h{p) 

1 -( v i-3)f 1 (p) + 1 -(3 + v 2 )f 2 (p) 
l J v 2 -3)h{p) + \{2,-v 2 )h{p) 



(4.23) 



Using (4.21), (4.22) these integrals can be done numerically without diffi- 



culties for any value of A. 
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As in the photon-photon case [29 , exact results can be obtained in the 
limits of small and large A. Using the known asymptotic properties of the 
function Gi(x) [501 EI] it is easy to show that, for small p, 



h(p) ~ ln(p) + | 7 + i In 3 + 0{p) 



Hp) 



(4.24) 



while for p — > oo 



MP) 

f2( P ) 



p 3 V 



(4.25) 



Using the leading terms of the expansions (4.25 ) in (4.23 ) yields the following 



results for the leading order terms of the amplitude ratios in the small A 
limit: 



a 32 
8^945 



A 2 + 0(A 4 ) , 



n^in(A) 



a 128 A 2 + 0(A 4 ) 



4vr 2835 



(4.26) 



Using the leading terms of the expansions (4.24) gives the asymptotic be- 
haviour for large A: 



n s %(A) = 


a 
8^ 


~MA) + 


4 2 , 
-7+ -ln3 - 
9 ' 3 


4, 

-In 2 + 
9 


! + o ( A-i) 


ft s l(A) = 


a 
8^ 


~MA) + 


4 2 , 


4, 

-In 2 + 
9 


| + 0(A-I) 


nSn(A) = 


a 

47T 


-gln(A) + 


8 7+ 4 ln3- 
9 ' 3 


8 ln2 + 
9 


| + o ( a-!, 


n^in(A) = 


a 

47T 


-gln(A) + 


8 4, 

-7+ -In 3 - 

9 ' 3 


8 ln2- 
9 


2 , 2 " 

27+°< A ~ S > 



(4.27) 
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Note that the leading logarithmic terms in (4.27) are independent of the po- 
larization choice. As in the case of the leading asymptotic growth for large 
field strength, eq. (4.6 ), their coefficients are directly related to the UV coun- 
terterms. This is another fact which is familiar from the photon-photon case. 
Fig. [3] shows the result of a num erical evaluation, using MATHEMATICA, 
of the parameter integral (4.23) for the amplitude ft® in (A). The small and 
large A approximations (4.26 ), (4.27) are also shown. 




Figure 3: Numerical plot of the amplitude ft® in (A) for A between and 20. 
The small and large A approximations are also shown. A global factor of ^ 
has been omitted. 



5 The electric field case 

The photon-graviton conversion process in a constant field seems to have 
been studied hitherto exclusively for the magnetic field case. The reason is, 
of course, the absence of evidence for the existence of very strong macro- 
scopic electric fields anywhere in nature. Still, in the near future it may be 
possible to reach field strengths of the order of E cr in the laboratory using 
optical or X-ray lasers |52| HkH 154"] . Thus it seems worthwhile to shortly dis- 
cuss also the electric case. Moreover, some of our results for the magnetic 
case carry over to the electric case simply by analytic continuation. Namely, 
electric-magnetic duality for this amplitude takes the following form (see 



eqs.(2.14) and (2.19)), 



21 



4®al,spin(^0,£) = n s ffi c i spin (<2>,£,0) 



B—*iE, uj^—iu) 
B^iE, u)— >— iu) 

(5.1) 



Despite of this formal duality, there is a structural difference between the 
electric and magnetic amplitudes. While the magnetic amplitude ratios are 
real for photon energies below the pair creation threshold, the integrands 
of the corresponding electric ones have poles in the Euclidean proper-time 
variable T = is, indicating the existence of an imaginary part for all photon 
energies. This again agrees with the photon-photon case |41| . and physically 
corresponds to the fact that both the photon/graviton and the electric field 
are capable of pair production, while in the magnetic case the field alone 
cannot induce an absorptive part. Although we are interested here only in 
the real parts, technically this implies that a calculation of the full electric 



amplitudes through numerical computation of the parameter integrals (3.5) 
is even more difficult than for the magnetic amplitudes. In particular, the 
rotation to Euclidean proper-time is much less useful here, even for u below 
threshold, since the rotated integrand has poles, both in the prefactor func- 
tions and in the exponent. We will therefore restrict our discussion to two 
limiting cases where we can directly draw on our results for the magnetic 
case, the zero photon energy case and the weak field case. 

5.1 Zero photon energy and arbitrary electric field strength 



For uj = 0, in eqs.(4.12) we have the magnetic amplitudes in an explicit 
form involving only the functions T, if), and Qh- Implementing the analytic 
continuation B — > iE is therefore straightforward. In particular, for large E 
the amplitude ratio shows the same logarithmic growth as we had found in 
the magnetic case (eq.(|4.6|)), 



U^ = 0,E) ™ --^—ln(E) , 



n&(a = o,£) -~in(£). 



127T 

a 
3ir 



(5.2) 
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5.2 Weak electric field and arbitrary photon energy 



Things become even simpler in the weak field limit. The approximation used 



in section 4.3 based on a weak field expansion and truncation to the leading 
order terms, removes the additional poles contained in the electric integrand, 
leading to a complete symmetry between the electric and magnetic cases. 
Since in this approxim ation the amplitude ratios depend only on the single 
variable A = %BCj, eqs.(5.1| turn into 



ft®" ( A 

scal,spin \ 



fr®-L f \ 

scal,spin \ 



-ECj 



-ECj 



ft®- 1 ( A 

scal,spin \ 



TT 11 A 

scal,spin \ 



-BCj 



-BCj 



(5.3) 



Therefore all the results of section 4.3 carry over to the electric case mutatis 
mutandis. 



6 Conclusions 

We have analyzed the parameter integral representations obtained in part 
1 for the one-loop photon-graviton amplitudes in a constant field at about 
the same level of detail as was previously achieved for the analogous QED 
photon-photon case [23 EH ED E2J [Ml HI]- In the purely magnetic case, 
we have shown that our representation is amenable to a direct numerical 
evaluation for photon/graviton energies below threshold, at arbitrary field 
strength. For weak magnetic fields, a one-parameter integral representation 
involving Airy functions has been obtained, and shown to be suitable to 
numerical evaluation. Closed-form results have been found for the zero 
energy limit. We have also transformed our results for the weak field and 
zero energy magnetic cases to corresponding results for the electric field case. 

The qualitative properties of the photon-graviton amplitudes turn out 
to be closely analogous to the ones of the corresponding photon-photon am- 
plitudes. In particular, in the magnetic case they have the same same pair 
creation thresholds, and a similar asymptotic growth for strong fields or 
large photon/graviton energies. This is, of course, not surprising, particu- 
larly considering that the photon-photon and photon-graviton amplitudes 
are connected by the gravitational Ward identity ((A. 13) of part 1). Apart 



from the pair creation thresholds (4.5), none of our results show any sub- 



stantial differences between the scalar and spinor loop cases. 
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From a quantitative point of view, our analysis of the one-loop photon- 
graviton amplitude can be summarized as follows. For small photon energies 
and large magnetic or electric field strengths, this ratio grows logarithmi- 
cally in the field strength (see (4.6) resp. (5.2)). For weak fields and large 
photon energies it grows logarithmically with the photon energy (see (4.27), 
(5.3)). In these limits it is clearly not possible to compensate the small 
prefactor a /4n for physically relevant values of the parameters. While we 
have not been able to perform a quantitative analysis in the whole two- 
parameter space of field strengths and energies, a compound effect of large 
field strength s and large energies appears to be exluded by the fact that the 
integrand of (3.5) contains u only in the combination lj 2 /B resp. Cj 2 /E. 
This is also borne out by the numerical results of section |4.1| Overall, we 
conclude that the magnitude of the one-loop contribution to this amplitude 
will not amount to more than a few percent of the tree level one for physically 
realistic values of the parameters. 

It should be mentioned, though, that there is also a qualitative difference 
to the tree-level amplitude. As has been stressed in |10| . the tree level 
photon-graviton conversion does, contrary to the photon-axion case, not 
lead to a dichroism effect for photon beams. This is because, according to 
(2.14), both photon polarization components have equal conversion rates. 
This symmetry does not extend to the one-loop level (except for the strong 
field limit). 

As we have seen, in the worldline formalism the calculation of the photon- 
graviton polarization tensor in a constant field is only moderately more 
difficult than the one of the photon-photon polarization tensor |38j. We 
expect that even the graviton-graviton case will be quite feasible. In a 
future sequel, we intend to analyze this case at the same level of the photon- 
graviton one, which then would make it possible to study the complete set 
of one-loop photon-graviton dispersion relations (1.4). 
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A Quadratic expansion of the Einstein-Maxwell 
theory 

The Einstein-Maxwell theory is described by 

%, A\ = j d D x ^ (±R - \f^F^ (A.l) 

where the metric g^ u has signature (— ,+,+,...,+), g = |det g^l, and k 2 = 
IGttGn. The spacetime dimension of interest is D = 4, but in this section we 
may as well keep it arbitrary. We expand g^ u = r/^ u +Kh^ u and = A^+a^. 
Then using the short-hand notation hn = d a h afl and h = rj^h^ one obtains 
the following quadratic approximation in the fluctuations (h^ u ,a^) around 
the background (77^, A^) 

% = Jd D xfyhrnh^-hnh + 2hd»hp + 2hPhr) + ^(aPnar + (dpa 



Kh llv (F^r a -- A t v F^f c 

(h 2 - \hl v )F 2 + h^h a ?F m F vP + {2h» a h\ - hh^)F^F/] } 



K 

T 



(A.2) 



In the second line of this expression we recognize the linear coupling %h iXV T llv 
of the graviton hnv with the stress tensor of the background electromagnetic 
field f^ v = W a F v a - \r)^ v F 2 . This tadpole vertex indicates that the non- 
trivial background stress tensor tends to curve the space. The third line in 
(A.2 ) gives instead the tree level graviton-photon mixing in the electromag- 
netic background. Using plane waves 

Ku{x) = e^e ikx , a a (x) = e a e ik * x (A.3) 

we get for this mixing term the vertex 

AS {2) = {2it) D 5 D {k + k 2 ) e^ u e a (-in) [F» a k u - (F ■ k)"r] ua + ~rT(F ■ k) 



(A.4) 

which appears in the path integrals as e % ^ s ^ (see eq. (1.2)). 
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The two-point functions (which we denote by II), in either coordinate or 
momentum space, are contained in S7 2 ) ( or m the quadratic part of the full 

effective action T( 2 ) = S( 2 ) + lo ° P ^ + • • •) as follows 



5 (2) 



d D x |i h^(x)n^ Xp (d)h Xp (x) + ^a a (x)IL a ^(d)a p (x) 

+h^{x)Ii^ a {d)a a {x) + | h^{x)f pv {x) 
d D k ( 1 1 



(2vr)^ \ 2 v ; Apv ; ' 2 

+V(^)n^' a (fc)a a (-A;) + ^V(fc)f^(-fc)} (A.5) 

where the Fourier transform of a field is given by h fll ,(x) = J i ^ e lkx h pu {k) , 

its inverse by h /iu (k) = J d D x e~ tkx h pv {x) and II(/c) = 11(9 — > — 

The equations of motion in term of these two-point functions then reacj^] 



Sa a (k) : U a ^(k)ap(-k) + U^ a (-k)h^(-k) = 

\f p ^-k) 

(A.6) 



*V(fc) : ^ u ' Xp (k)h Xp (-k) + U^ a (k)a a (-k) = --f^(-k) 



and, in particular, one obtains from (A. 2) 



n-tr'ee(^) ~~ k a kP k 2 T] a ^ 



iCT(fc) = -y C^{k) (A.7) 



with C pu ' a (k) as in (1.2). Also, one gets 

„2 



-^tree 



(fc) = ^- (VV P - ^«V P + ^V*)) - l -(t V k X k p + k p k u V xp ) 

IIS. 

(A.8) 



i f 5 i - 

i 



+ -(r] px k u k p + r] uX k p k p + rj pp k u k x + r] vp k^k x ) + F terms. 



Finally, note that a constant T pu gives T^(A;) = (2ir) 2 5 D (k)T pu . 

3 A minus sign in the euclidean formulas (A. 11) of paper I is correctly taken into account 
by the Wick rotation, as one can easily check by looking at the euclidean tree level Maxwell 
action S — J d D x jF^F^", which indeed is positive definite in euclidean space (as it 
should) . 
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